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Abstract: The lepton mixing angle ^13, the only unknown angle in the standard three- 
flavor neutrino mixing scheme, is finally measured by the recent reactor and accelerator 
neutrino experiments. We perform a combined analysis of the data coming from T2K, 
MINOS, Double Chooz, Daya Bay and RENO experiments and find sin^ 26*13 = 0.096 ± 
0.013(ib0.040) at 1 cj (3 a) CL and that the hypothesis ^13 = is now rejected at a 
significance level of 7.7 a. We also discuss the near future expectation on the precision of 
the ^13 determination by using expected data from these ongoing experiments. 
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1 Introduction 

The accelerator search for fg appearance [1]^ and the precision measurement of reactor 
neutrino disappearance [2, 3] are both viable ways to measure ^13, which has been, until 
very recently, the unique unknown mixing angle of the lepton flavor mixing matrix [4]. 
It must be stressed that the experimental redundancy for measuring 6*13 may be justified 
because of the complementary nature of the two types of experiments, as discussed, for 
example, in [3, 5]. While the reactor experiments provide a clean measurement of ^13 which 
is free from degeneracy [6-8], the accelerator measurement can enjoy the interplay with 
the CP phase 5cf, which connotes the possibility of extension of the experiment to an 
upgraded phase to search for CP violation. 

It is very fortunate to see that the era of simultaneous measurement of ^13 by accel- 
erator and reactor has just arrived. In June of 2011 the T2K group reported six clean 
events of Vg appearance, implying 2.5 a indication for non-zero ^13 [9] with a best fit value 
comparable to the CHOOZ limit [10] (see also [11-13]). It was soon followed by the MINOS 
collaboration which reported also indication of non-zero ^13 [14]. At the end of 2011, one 
of the reactor ^13 experiments. Double Chooz [15], reported their first result, constraining 
(9i3 to a range sin^ 2(9i3 = 0.086 ±0.051 at 68 % CL [16, 17]. ^ 

Very recently we were surprised by the announcement of another two reactor ^13 exper- 
iments. First, Daya Bay [18], which reported a measurement of ^13 as accurate as sin^ 2^13 
= 0.092 ±0.017 at 68 % CL. The significance level for non-zero ^13 obtained by them is 5.2 

^For an updated version see http://neutrino.kek.jp/jhfnu/loi/loi.v2.030528.pdf 
^For the official Double Chooz results see http://doublechooz.in2p3.fr/Status_and_News/status_ 
and_news .php. 
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a [19]. Second, RENO[20, 21], which reported sin^ 2(9i3 = 0.113 ±0.023 at 68 % CL, ex- 
cluding a non-zero ^13 at 4.9 a [22]. Though still limited by both statistics and systematics 
(except for Daya Bay whose systematics is already small), these results, together with the 
aforementioned accelerator data, constitutes the most valuable information on ^13 to date. 
Therefore, we believe that it is a meaningful step to attempt a combined analysis of these 
data set. 

The issue of possible non-zero ^13 has been discussed in the context of global analyses 
which include the solar and the atmospheric neutrino data even before [23-25] or after [26, 
27] the T2K result [9]. However, given the current precision on the determination of ^13, 
such global fits seem unnecessary in this specific context. Hence, in this paper we restrict 
ourselves to a combined analysis of the accelerator and the reactor ^13 experiments only. 

2 Analysis details 

We analyze the available accelerator data from T2K [9] and MINOS [14] in the — )• Ue 
appearance channel in combination with the very recent Double Chooz [16, 17], Daya 
Bay [19] and RENO [22] reactor data in the Ve z^e disappearance channel. We will also 
make some prognostication to the near future. The simulations were performed using a 
modified version of GLoBES [28]. 

2.1 Accelerator experiments: T2K and MINOS 

The T2K experiment uses a narrow 2.5° off-axis beam generated at J-PARC in Tokai 
which is directed to the Super-Kamiokande detector of fiducial mass 22.5 kt located in 
Kamioka 295 km away from J-PARC. In order to reproduce the T2K allowed region in the 
sin^ 2^13 - 5cp plane, reported in figure 6 of ref. [9], we have simulated the T2K signal in 
the Vfj_ — >• Ve appearance channel in a similar way as done in ref. [5]. We took the neutrino 
fluxes from the letter of intent of the Hyper-Kamiokande project [29] and the background 
from [9]. The cross sections and energy dependent efficiencies for charged current quasi- 
elastic (CC-QE) and non quasi-elastic (CC-NQE) events are simulated in a similar manner 
as in [5] to reproduce the energy spectra given in [29]. 

Energy smearing and the consequent migration of events were taken into account in 
our calculations by using a Gaussian energy resolution function with width 85 (130) MeV 
for CC-QE (CC-NQE) events. For CC-NQE events, following the procedure described in 
the Appendix of [5], a shift of 350 MeV was introduced in the Gaussian smearing function 
in order to take into account the significant difference between true and reconstructed 
neutrino energy. In reproducing the current T2K result we assumed 1.43 x 10^° POT and 
23% systematic uncertainly in the absolute normalization. 

The MINOS experiment uses the NuMI beamline and operates with a near detector 
located on-site at Fermilab, and a far detector located 735 km away in the Soudan Under- 
ground Laboratory. The near (far) detector consists of 0.98 kt (5.4 kt) of alternating layers 
of steel and plastic scintillator. In order to reproduce the MINOS allowed region in the 
sin^ 2^13 — 5cp plane, given in figure 3 of ref. [14], we have simulated the Ve signal using 
the same procedure as in ref. [30] but with the background and systematic uncertainties 
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taken from [14]. We assumed a total exposure of 8.2 X 10^0 POT, but a tuning of the 
normalization was needed in order to obtain the correct number of signal events. 

2.2 Reactor experiments: Double Chooz, Daya Bay and RENO 

Double Chooz (DC) is a reactor antineutrino oscillation experiment [15] based on the 
CHOOZ-B Nuclear Power Station. The experiment is a double detector apparatus (each 
detector with a fiducial volume of 10.3 m^) based on liquid scintillator, though until 2013 
they will be taking data only with their far detector located at 1.05 km from the two 4.27 
GWth reactor cores. 

To simulate the Ve Vf. disappearance reported by DC collaboration in refs. [16, 17] 
we have performed a calculation based on the far detector specification and reactor fuel 
composition given in ref. [15], with systematic uncertainties, background and efficiency, 
and other additional information according to [16, 17]. 

Before analyzing the real data, we first tried to reproduce the expected visible energy 
spectra obtained by the Monte Carlo (MC) simulations of the DC collaboration, in the 
absence and presence of oscillation shown (respectively, by the blue dotted and red solid 
histograms) in figure 3 of [16]. Indeed, in our attempt to reproduce the visible energy 
spectra, we have noticed that these spectra exhibit significant distortions if compared to 
the corresponding spectra as a function of the true prompt energy, which, of course, can 
not be measured directly. 

In order to mimic such a rather strong distortions, which are due to various effects 
taken into account in the MC simulations by the DC collaboration, we first introduce an 
energy smearing effects using a Gaussian energy resolution function with a width cfe = 
l2%^{E/MeV) + 0.15 MeV. We note that due to the 2nd term in erg, we can reproduce 
rather well the spectra after taking into account the additional corrections described below. 
We, however, stress that the inclusion or omission of the 2nd term in a^; does not alter 
much the allowed parameter region of sin'^ 2^13 and (5cp presented in this paper, though it 
affects the Xmin values. 

In addition to the energy smearing we have further taken into account, in an approx- 
imate way, two kinds of corrections which were actually introduced in the analysis by the 
DC collaboration [31]^ in order to understand their data. The first one is a non linearity 
correction. This is based on the energy calibration by using several sources performed 
by the DC collaboration. Roughly speaking, the observed visible energies (or to be more 
precise the number of photoelectrons) tend to be overestimated (underestimated) for en- 
ergy larger (smaller) than ~ 1.5 MeV for up to a few percent, when compared to the ones 
predicted by MC simulations. Note that the correction is energy dependent, see [31]. The 
second correction is the one based on the Z-dependence calibration, which shows that the 
observed energy tends to be underestimated when the neutrino event occurs in the region 
far from the center of the detector, for up to a few percent [31]. We note that, after taking 
into account these two corrections in addition to the energy smearing, we can reproduce 
reasonably well the energy spectra shown in figure 3 of [16]. 

^Relevant information to be available at http://doublechooz.in2p3.fr/Status_and_News/status_ 
and_news .php. 
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Daya Bay (DB) experiment measures Ve from six 2.9 GWth reactors grouped into three 
pairs of nuclear power plants (NPP) using six detectors deployed in a near-far arrangement 
allowing to compare rates at various baselines. There are two near detectors at 364 m from 
the Daya Bay NPP, one near detector at 480 m (528 m) from Ling Ao (Ling Ao-II) NPP and 
three far detectors at 1912 m (1540 m) from Daya Bay (Ling Ao and Ling Ao-II) NPP [19]. 
Each one of these identical detectors is made of a 5 m diameter cylindrical stainless steel 
vessel which holds a 3.1 m diameter inner vessel and a 4 m diameter outer vessel. The inner 
vessel holds 20 ton of Gadolinium-doped liquid scintillator (the target) which is shielded 
from the 20 ton liquid scintillator in the outer acrylic vessel by 37 ton of mineral oil. We 
have implemented in our simulation the systematic uncertainties and the background in 
accordance with [19]. In this paper, for DB, we do not consider the energy spectrum but 
restrict ourselves to the rates only analysis, as done in ref. [18], even for our near future 
analysis to be discussed in section 4. The reason is that DB's systematic uncertainty for 
the rate only analysis is already quite small and their data are still statistically limited and 
hence it seems that the spectrum information would not play an important role for the 
time being. 

RENO is the reactor experiment which receives neutrinos from the YongGwang Nuclear 
Power Plant located 400 km from Seoul in which six 2.8 GWth reactors are lined up. They 
use two 16 t liquid scintillator identical detectors, the near (far) detector located at roughtly 
300 m (1.3 km) from the reactors. RENO has been taking data with both detectors since 
August 2011. In order to simulate RENO v^, disappearance signal we performed rate-only 
analysis, using the background, energy resolution and the systematic uncertainties given 
in ref. [21, 22]. 

In DC, DB and RENO simulations we have used the new reactor antineutrino flux 
calculations [32, 33]. This has little impact on our results, since DC is normalized to 
Bugey-4 cross section and the other reactor experiments have near detectors. 

3 Analysis results: current status 

3.1 Combining accelerator and reactor data 

Before combining the accelerator and reactor neutrino data we have verified that we are 
able to reproduce quite well the individual result of each experiment T2K [9], MINOS [14], 
DC [17], DB [19] and RENO [22]. Here we present our combined analysis of these experi- 
ments. 

In figure 1 we show the allowed region obtained in our combined analysis. The yellow, 
orange, and red bands correspond, respectively, to 68%, 95%, 99% CL regions for 2 degree 
of freedom (dof). We also show the behavior of Ax^ for 1 dof as a function of sin^ 2^13 
(attached to top), and as a function of 5cp (attached to right side) in each panel. For T2K, 
DB and RENO, only the total rate was considered, whereas for MINOS (DC) we used data 
of 7 (18) energy bins. However, we have checked that T2K allowed region does not change 
much if we also take into account the spectrum information. In our fit we have explicitly 
assumed one of the mass hierarchies (normal or inverted) as input and varied sin^ 2^23 
and I Am32 1 , imposing Gaussian priors based on the atmospheric neutrino experiments [34] 
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Figure 1. Allowed region in sin^ 26I13 - 6cp plane for T2K, MINOS, Double Chooz (DC), Daya 
Bay (DB) and RENO combined at 68%, 95 % and 99% CL for 2 dof, assuming normal (left panel) 
or inverted (right panel) mass hierarchy. We also show the behavior as a function of sin'^ 2^13 
(top) and as a function of Sqp (right) in each case. As a reference we also show the 90%CL exclusion 
limit from CHOOZ [10]. 

and MINOS [35] results. We observe that if we combine only T2K and DC (not shown in 
figure 1), our allowed regions agree very well with the result shown in [17] for the same 
fixed values of sin^ 2^23 and | Am^g | • 

We conclude that at 95% CL, the allowed range for ^13 is given as 0.070 < sin^ 2^13 < 
0.122 irrespectively of the mass hierarchy for 1 dof. In the case of normal (inverted) 
mass hierarchy, the best fit point is given by sin^ 2^13 = 0.096 (sin^ 2^13 = 0.096) and 
6cp = 0.977r {6cp = -0.147r) which correspond to Xmin/(24 - 2) = 1.57 (1.55). At the 
moment, there is not much significance in the preferred value of 6cp- We can also see the 
contribution of each individual experiment to the determination of sin^ 2^13. Currently 
DB is the most powerful experiment in constraining sin^ 2^13 from both ends. Before DB 
and RENO announced their results T2K was the most effective experiment in excluding 
a vanishing value of sin^ 2^13 , but allowed for higher values of sin^ 2^13 than MINOS and 
DC. The combination of the five experiments can now exclude sin'^ 2^13 = at 7.7 a CL, 
irrespectively of the mass hierarchy. 

3.2 Potential hint on CP violation 

It was proposed in [36] that hints of CP violation could be obtained by combining accel- 
erator and reactor measurements. In this method, determining sgn(sin(5cp) is essentially 
the goal to reach. At this moment, however, change in Ax^ is quite mild, < 1.0, as 6cp 
is varied, as we can see from figure 1. Clearly, it is not possible to make any definitive 
statements about which sign of sinJcP is preferred. Nevertheless, we may say that the 
region sgn(sin(5cp) > (sgn(sin(5cp) < 0) is slightly preferred in the normal (inverted) 
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mass hierarchy case. This tendency could become clearer by future accumulation of the 
data, as shown in section 4. 

We note that another hint of sgn(sin(5cp) is provided by the three flavor analysis of 
the SK atmospheric neutrino data [37, 38], which indicates, though mildly, the negative 
sin^cp region, which is more prominent in the inverted hierarchy case. We believe that 
the issue of preferred region of 5cp deserves careful watching with accumulation of various 
experimental data in the future. 

4 Expectation: one year from now 

We now make some predictions for the possible situation of ^13 in the near future, about 
one year from now. For definiteness, in our simulations for the future expectation, we 
assume the true parameters to be our best fit value sin^ 2^13 = 0.096, 6cp = O.OTvr and 
the normal hierarchy scheme, though we confirmed that the results do not change much 
even if the inverted hierarchy (with the corresponding best fit values) was assumed. We do 
not consider MINOS in our predictions because the impact of the improvement of MINOS 
sensitivity to ^13 appears to be limited. We include the energy spectrum information in the 
analysis of T2K. We used the same priors as before for |Am|2| and sin^2023- While this 
may seem too conservative, these uncertainties mainly have an effect on the upper bound 
on sin^ 2^13. 

We take the same systematic uncertainties and the backgrounds claimed by the ex- 
periments, as in the previous section. (For T2K and DC, we also considered the case with 
reduced systematic uncertainties, see the footnote 5 and text below.) We assume that 
DC have been taking data since April 2011 with averaged 77.5% data taking efficiency 
for physics and 76% reactor power efficiency to take into account reactor off periods. For 
RENO, we set the data taking to start at August 2011 and used the efficiencies and DAQ 
live times (proportionally) given in ref. [20]. 

We assume T2K will resume its operation in January 2012 with their proposed inte- 
grated luminosity of 10^^ POT/year.'^ For concreteness, we assume the present configu- 
ration of DB throughout the year. As they increase the number of detectors in the far 
experimental hall one can shift our results by the appropriate number of days/months to 
roughly account for that. 

On the left panel of figure 2 we show the expected 1 a uncertainty on the determination 
of sin^ 2^13 as a function of time, for the different experiments. We employ the following 
color code for the bands: pink for DC, green for RENO, light blue for T2K, blue for DB 
and yellow for the combination. 

We observe that, at this moment, DB with six detectors is the most powerful experi- 
ment among the four and dominates the final combined result. RENO, the next powerful 
one, would have dominated the combination at the very beginning of this year. However, 
as soon as their result becomes dominated by systematic uncertainties they cannot improve 

* We know that this assumption no longer holds. However, we remain with it because we do not know 
for sure the real situation of the T2K experiment in 2012. Therefore, as far as T2K is concerned, our 
predictions can be viewed as optimistic. 
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Figure 2. On the left panel, we show the expected 1 a uncertainty on sin^ 2^13 for the case 
where the true value of sin^ 2^13 = 0.096 and 5 = 0.977r (current best fit for the normal hierarchy) 
as a function of the months in 2012 for DC, RENO, T2K, DB as well as the combined case. On 
the right panel, we show the expected 1-3 a uncertainties on sin^ 2^13 as a function of time for the 
same input but only for the combined case. On the left panel we indicate the effect of the improved 
systematics considered for T2K and DC by the dashed and dash-dotted lines, respectively, whereas 
their impact on the combined analysis is indicated by the dotted line on the left panel and by the 
dotted, dashed and solid lines on the right panel (see legends in the plots). In fitting, the hierarchy 
is assumed to be unknown. 



their sensitivity much. To do that they will have to improve their systematics and/or do 
a spectrum analysis. DB, due to its high reactor power, overall detector mass and smaller 
systematic uncertainty, quickly becomes the most powerful experiment and remains so 
throughout the year. DC sensitivity can not improve much with only the far detector. 
Regarding the bound on sin^ 2^13 from below in 2012, T2K is comparable to DC, however 
both will not reach the discriminability of RENO and DB. 

We note that although a reduction of systematic uncertainties by ~ 30-50% for T2K 
and DC ^ have some impact on the individual results of these experiments, it does not 
essentially affect the accuracy of determination of ^13 based on the combined data in the 
present circumstances. On the other hand, a possible reduction of RENO backgrounds can 
fairly affect the final sensitivity to ^13. See figure 2. The improvement in T2K is less visible 
probably because it is still dominated by statistical error. 

On the right panel of figure 2, we show the la — 2>a uncertainty regions for the deter- 
mination of sin^ 2^13 as a function of time for all experiments combined. Here the yellow, 

^ For the DC experiment we use for each systematic uncertainty the best value between the one quoted 
in their proposal and the one presented in their latest paper [16]. We also apply a 30% reduction of the 
number of their background events. For T2K we assume about 50% reduction of systematic uncertainties, 
arbitrarily re-scaling the current normalization uncertainty to 10%. For RENO, we only contemplate the 
possibility of reducing the backgrounds by 40% at the cost of decreasing the near and far detector signal 
in 10% and 3%, respectively, while for DB we do not consider any reduction of the systematics due to the 
already extraordinarily low systematics of this experiment. 
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orange and red bands correspond, respectively, to la, 2a, and 3a regions. From this anal- 
ysis we conclude that within 1 year, the uncertainty on the determination of sin^ 2^13 at 
la CL may be reduced from 0.013 to ~ 0.005. We have verified that this is also true for 
the case where the true mass hierarchy is the inverted one. 

At the same time, the hypothesis of a vanishing sin^ 2^13 could be rejected at a level 
of very high significance. We have verified, under the above stated assumptions, that by 
the middle (end) of 2012, the sin^ 2&13 = hypothesis could be rejected with a significance 
larger than 11 a (14 a), if the future data is consistent with the current best fit point. 
We confirm that the impact of the improvement of systematics for T2K and DC on the 
combined analysis is quite small also for the 2 and 3 a regions (see dashed and solid lines 
on the right panel of figure 2). 

In figures 3 and 4, using the same format as in figure 1, we show the allowed region 
in the sin^ 2^13 — 6cp plane expected in June and in December of 2012, respectively, that 
could be achieved by combining T2K, MINOS, DC, DB and RENO data. As input, we used 
the best fit point for the normal hierarchy and fitted for each normal and inverted mass 
hierarchy. As expected, from the right panel of figure 2, the impact of the reduction of the 
systematic uncertainties we considered in this work on the determination of the parameter 
regions, as well as in the behavior of Ax^ (indicated by the dashed lines in figures 3 and 
4), is quite small as far as the results expected in the near future (~ 1 year) are concerned. 

Finally, we note that at the end of this year the combined Ax^ for different values 
of 5cp is expected to be ~ 1-4, depending on the fitted hierarchy. This might be used 
as a hint on which region of 6cp is preferred, but this still will not be strong enough to 
definitively pin down the value of 5cp with high significance. For future prospects on the 
reactor-accelerator combined method, see also [39]. 

5 Conclusion 

We performed a combined analysis of the currently available accelerator and reactor data 
which provide a very significant evidence of non-zero ^13. Being outside of the experi- 
mental collaborations our simulation may be incomplete by lack of sufficient information 
on backgrounds, systematic uncertainties, efficiencies, etc. However, we believe that we 
did a reasonable job and our results serve as an independent confirmation of the analyses 
provided by the experimental groups. 

It is encouraging to see that the confidence level for non-zero sin^ 2^13 now reaches 
~ 7.7 a thanks, in particular, to DB and RENO, in addition to T2K, MINOS and DC 
experiments. Still in this year, we will have indisputable evidence for non zero ^13. We 
predict that if the future data continues to be compatible with the current best fit value 
sin^ 2^13 = 0.096, by the middle (end) of 2012 sin^ 26*13 will be known within ±0.007 
(±0.005) at 68% CL. Finally, we also studied the impact of the possible reduction of the 
systematic uncertainties for DC (by 30%) and T2K (roughly by half), as well as for RENO 
backgrounds (by 40%) while forfeiting signal efficiency (3% in the far detector). We have, 
however, found that the reduction of errors do not affect in a significant way the combined 
sensitivity at the end of 2012. 
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Figure 3. Predicted allowed region in the sin^ 26*13 - (^cp plane for T2K, MINOS, DC, DB and 
RENO combined at 68%, 95 % and 99% CL for 2 dof in the middle (June) of 2012, assuming normal 
(left panel) or inverted (right panel) mass hierarchy and as input the normal hierarchy best fit point 
of our current analysis. 
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